Vesicular trafficking is a key determinant of the statin response in acute myeloid leukemia.

Cholesterol homeostasis has been proposed as one mechanism contributing to chemoresistance in AML and hence, inclusion of statins in therapeutic regimens as part of clinical trials in AML has shown encouraging results. Chemical screening of primary human AML specimens by our group led to the identification of lipophilic statins as potent inhibitors of AMLs from a wide range of cytogenetic groups. Genetic screening to identify modulators of the statin response uncovered the role of protein geranylgeranylation and of RAB proteins, coordinating various aspect of vesicular trafficking, in mediating the effects of statins on AML cell viability. We further show that statins can inhibit vesicle-mediated transport in primary human specimens, and that statins sensitive samples show expression signatures reminiscent of enhanced vesicular trafficking. Overall, this study sheds light into the mechanism of action of statins in AML and identifies a novel vulnerability for cytogenetically diverse AML.

High frequency of germline RUNX1 mutations in patients with RUNX1-mutated AML.

RUNX1 is mutated in ∼10% of adult acute myeloid leukemia (AML). Although most RUNX1 mutations in this disease are believed to be acquired, they can also be germline. Indeed, germline RUNX1 mutations result in the well-described autosomal-dominant familial platelet disorder with predisposition to hematologic malignancies (RUNX1-FPD, FPD/AML, FPDMM); ∼44% of affected individuals progress to AML or myelodysplastic syndromes. Using the Leucegene RUNX1 AML patient group, we sought to investigate the proportion of germline vs acquired RUNX1 mutations in this cohort. Our results showed that 30% of RUNX1 mutations in our AML cohort are germline. Molecular profiling revealed higher frequencies of NRAS mutations and other mutations known to activate various signaling pathways in these patients with RUNX1 germline-mutated AML. Moreover, 2 patients (mother and son) had co-occurrence of RUNX1 and CEBPA germline mutations, with variable AML disease onset at 59 and 27 years, respectively. Together, these data suggest a higher than anticipated frequency of germline RUNX1 mutations in the Leucegene cohort and further highlight the importance of testing for RUNX1 mutations in instances in which allogeneic stem cell transplantation using a related donor is envisioned.

Genetic characterization of ABT-199 sensitivity in human AML.

Acute myeloid leukemias (AML) with mutations in the NPM1 gene (NPM1c+) represent a large AML subgroup with varying response to conventional treatment, highlighting the need to develop targeted therapeutic strategies for this disease. We screened a library of clinical drugs on a cohort of primary human AML specimens and identified the BCL2 inhibitor ABT-199 as a selective agent against NPM1c+ AML. Mutational analysis of ABT-199-sensitive and -resistant specimens identified mutations in NPM1, RAD21, and IDH1/IDH2 as predictors of ABT-199 sensitivity. Comparative transcriptome analysis further uncovered BCL2A1 as a potential mediator of ABT-199 resistance in AML. In line with our observation that RAD21 mutation confers sensitivity to ABT-199, we provide functional evidence that reducing RAD21 levels can sensitize AML cells to BCL2 inhibition. Moreover, we demonstrate that ABT-199 is able to produce selective anti-AML activity in vivo toward AML with mutations associated with compound sensitivity in PDX models. Overall, this study delineates the contribution of several genetic events to the response to ABT-199 and provides a rationale for the development of targeted therapies for NPM1c+ AML.

Mubritinib Targets the Electron Transport Chain Complex I and Reveals the Landscape of OXPHOS Dependency in Acute Myeloid Leukemia.

To identify therapeutic targets in acute myeloid leukemia (AML), we chemically interrogated 200 sequenced primary specimens. Mubritinib, a known ERBB2 inhibitor, elicited strong anti-leukemic effects in vitro and in vivo. In the context of AML, mubritinib functions through ubiquinone-dependent inhibition of electron transport chain (ETC) complex I activity. Resistance to mubritinib characterized normal CD34+ hematopoietic cells and chemotherapy-sensitive AMLs, which displayed transcriptomic hallmarks of hypoxia. Conversely, sensitivity correlated with mitochondrial function-related gene expression levels and characterized a large subset of chemotherapy-resistant AMLs with oxidative phosphorylation (OXPHOS) hyperactivity. Altogether, our work thus identifies an ETC complex I inhibitor and reveals the genetic landscape of OXPHOS dependency in AML.

Chemogenomic Landscape of RUNX1-mutated AML Reveals Importance of RUNX1 Allele Dosage in Genetics and Glucocorticoid Sensitivity.

Purpose:RUNX1-mutated (RUNX1mut) acute myeloid leukemia (AML) is associated with adverse outcome, highlighting the urgent need for a better genetic characterization of this AML subgroup and for the design of efficient therapeutic strategies for this disease. Toward this goal, we further dissected the mutational spectrum and gene expression profile of RUNX1mut AML and correlated these results to drug sensitivity to identify novel compounds targeting this AML subgroup.Experimental Design: RNA-sequencing of 47 RUNX1mut primary AML specimens was performed and sequencing results were compared to those of RUNX1 wild-type samples. Chemical screens were also conducted using RUNX1mut specimens to identify compounds selectively affecting the viability of RUNX1mut AML.Results: We show that samples with no remaining RUNX1 wild-type allele are clinically and genetically distinct and display a more homogeneous gene expression profile. Chemical screening revealed that most RUNX1mut specimens are sensitive to glucocorticoids (GCs) and we confirmed that GCs inhibit AML cell proliferation through their interaction with the glucocorticoid receptor (GR). We observed that specimens harboring RUNX1 mutations expected to result in low residual RUNX1 activity are most sensitive to GCs, and that coassociating mutations as well as GR levels contribute to GC sensitivity. Accordingly, acquired glucocorticoid sensitivity was achieved by negatively regulating RUNX1 expression in human AML cells.Conclusions: Our findings show the profound impact of RUNX1 allele dosage on gene expression profile and glucocorticoid sensitivity in AML, thereby opening opportunities for preclinical testing which may lead to drug repurposing and improved disease characterization. Clin Cancer Res; 23(22); 6969-81. ©2017 AACR.

UBAP2L is amplified in a large subset of human lung adenocarcinoma and is critical for epithelial lung cell identity and tumor metastasis.

The ubiquitin-associated protein 2-like (UBAP2L) gene remains poorly studied in human and mouse development. UBAP2L interacts with the Polycomb group protein B lymphoma Mo-MLV insertion region 1 homolog (BMI1) and determines the activity of mouse hematopoietic stem cells in vivo Here we show that loss of Ubap2l leads to disorganized respiratory epithelium of mutant neonates, which die of respiratory failure. We also show that UBAP2L overexpression leads to epithelial-mesenchymal transition-like phenotype in a non-small cell lung carcinoma (NSCLC) cell line. UBAP2L is amplified in 15% of human primary lung adenocarcinoma specimens. Such patients express higher levels of UBAP2L and show a reduction in survival when compared with those who do not have this gene amplification. Supporting a possible role for UBAP2L in lung tumor progression, NSCLC cells engineered to express low levels of this gene produce much smaller tumors in vivo than wild-type control cells. Together, these results suggest that UBAP2L contributes to epithelial lung cell identity in mice and that it plays an important role in human lung adenocarcinoma.-Aucagne, R., Girard, S., Mayotte, N., Lehnertz, B., Lopes-Paciencia, S., Gendron, P., Boucher, G., Chagraoui, J., Sauvageau, G. UBAP2L is amplified in a large subset of human lung adenocarcinoma and is critical for epithelial lung cell identity and tumor metastasis.

MiSTIC, an integrated platform for the analysis of heterogeneity in large tumour transcriptome datasets.

Genome-wide transcriptome profiling has enabled non-supervised classification of tumours, revealing different sub-groups characterized by specific gene expression features. However, the biological significance of these subtypes remains for the most part unclear. We describe herein an interactive platform, Minimum Spanning Trees Inferred Clustering (MiSTIC), that integrates the direct visualization and comparison of the gene correlation structure between datasets, the analysis of the molecular causes underlying co-variations in gene expression in cancer samples, and the clinical annotation of tumour sets defined by the combined expression of selected biomarkers. We have used MiSTIC to highlight the roles of specific transcription factors in breast cancer subtype specification, to compare the aspects of tumour heterogeneity targeted by different prognostic signatures, and to highlight biomarker interactions in AML. A version of MiSTIC preloaded with datasets described herein can be accessed through a public web server (http://mistic.iric.ca); in addition, the MiSTIC software package can be obtained (github.com/iric-soft/MiSTIC) for local use with personalized datasets.

Expression of immunoproteasome genes is regulated by cell-intrinsic and -extrinsic factors in human cancers.

Based on transcriptomic analyses of thousands of samples from The Cancer Genome Atlas, we report that expression of constitutive proteasome (CP) genes (PSMB5, PSMB6, PSMB7) and immunoproteasome (IP) genes (PSMB8, PSMB9, PSMB10) is increased in most cancer types. In breast cancer, expression of IP genes was determined by the abundance of tumor infiltrating lymphocytes and high expression of IP genes was associated with longer survival. In contrast, IP upregulation in acute myeloid leukemia (AML) was a cell-intrinsic feature that was not associated with longer survival. Expression of IP genes in AML was IFN-independent, correlated with the methylation status of IP genes, and was particularly high in AML with an M5 phenotype and/or MLL rearrangement. Notably, PSMB8 inhibition led to accumulation of polyubiquitinated proteins and cell death in IPhigh but not IPlow AML cells. Co-clustering analysis revealed that genes correlated with IP subunits in non-M5 AMLs were primarily implicated in immune processes. However, in M5 AML, IP genes were primarily co-regulated with genes involved in cell metabolism and proliferation, mitochondrial activity and stress responses. We conclude that M5 AML cells can upregulate IP genes in a cell-intrinsic manner in order to resist cell stress.

Chemo-genomic interrogation of CEBPA mutated AML reveals recurrent CSF3R mutations and subgroup sensitivity to JAK inhibitors.

In this study, we analyzed RNA-sequencing data of 14 samples characterized by biallelic CEBPA (CEBPA(bi)) mutations included in the Leucegene collection of 415 primary acute myeloid leukemia (AML) specimens, and describe for the first time high frequency recurrent mutations in the granulocyte colony-stimulating factor receptor gene CSF3R, which signals through JAK-STAT proteins. Chemical interrogation of these primary human specimens revealed a uniform and specific sensitivity to all JAK inhibitors tested irrespective of their CSF3R mutation status, indicating a general sensitization of JAK-STAT signaling in this leukemia subset. Altogether, these results identified the co-occurrence of mutations in CSF3R and CEBPA in a well-defined AML subset, which uniformly responds to JAK inhibitors and paves the way to personalized clinical trials for this disease.

GPR56 identifies primary human acute myeloid leukemia cells with high repopulating potential in vivo.

Acute myeloid leukemia (AML) is a genetically heterogeneous hematologic malignancy, which is initiated and driven by a rare fraction of leukemia stem cells (LSCs). Despite the difficulties of identifying a common LSC phenotype, there is increasing evidence that high expression of stem cell gene signatures is associated with poor clinical outcome. Identification of functionally distinct subpopulations in this disease is therefore crucial to dissecting the molecular machinery underlying LSC self-renewal. Here, we combined next-generation sequencing technology with in vivo assessment of LSC frequencies and identified the adhesion G protein-coupled receptor 56 (GPR56) as a novel and stable marker for human LSCs for the majority of AML samples. High GPR56 expression was significantly associated with high-risk genetic subgroups and poor outcome. Analysis of GPR56 in combination with CD34 expression revealed engraftment potential of GPR56(+)cells in both the CD34(-)and CD34(+)fractions, thus defining a novel LSC compartment independent of the CD34(+)CD38(-)LSC phenotype.

The transcriptomic landscape and directed chemical interrogation of MLL-rearranged acute myeloid leukemias.

Using next-generation sequencing of primary acute myeloid leukemia (AML) specimens, we identified to our knowledge the first unifying genetic network common to the two subgroups of KMT2A (MLL)-rearranged leukemia, namely having MLL fusions or partial tandem duplications. Within this network, we experimentally confirmed upregulation of the gene with the most subtype-specific increase in expression, LOC100289656, and identified cryptic MLL fusions, including a new MLL-ENAH fusion. We also identified a subset of MLL fusion specimens carrying mutations in SPI1 accompanied by inactivation of its transcriptional network, as well as frequent RAS pathway mutations, which sensitized the leukemias to synthetic lethal interactions between MEK and receptor tyrosine kinase inhibitors. This transcriptomics-based characterization and chemical interrogation of human MLL-rearranged AML was a valuable approach for identifying complementary features that define this disease.

Gene expression profile in long-term non progressor HIV infected patients: in search of potential resistance factors.

Long-term non-progressors (LTNP) represent a minority (1-5%) of HIV-infected individuals characterized by documented infection for more than 7-10 years, a stable CD4+ T cell count over 500/mm(3) and low viremia in the absence of antiretroviral treatment. Protective factors described so far such as the CCR5delta32 deletion, protective HLA alleles, or defective viruses fail to fully explain the partial protection phenotype. The existence of additional host resistance mechanisms in LTNP patients was investigated here using a whole human genome microarray study comparing gene expression profiles of unstimulated peripheral blood mononuclear cells from LTNP patients, HIV-1 infected patients under antiretroviral therapy with CD4+ T cell levels above 500/mm(3) (ST), as well as healthy individuals. Genes that were up- or downregulated exclusively in LTNP, ST or in both groups in comparison to controls were identified and classified in functional categories using Ingenuity Pathway Analysis. ST and LTNP patient groups revealed distinct genetic profiles, regarding gene number in each category and up- or downregulation of specific genes, which could have a bearing on the outcome of each group. We selected some relevant genes to validate the differential expression using quantitative real-time qRT-PCR. Among others, we found several genes related to the canonical Wnt/beta-catenin signaling pathway. Our results identify new possible host genes and molecules that could be involved in the mechanisms leading to the slower progression to AIDS and sustained CD4+ T cell counts that is peculiar to LTNP patients.

Impact of genomic polymorphisms on the repertoire of human MHC class I-associated peptides.

For decades, the global impact of genomic polymorphisms on the repertoire of peptides presented by major histocompatibility complex (MHC) has remained a matter of speculation. Here we present a novel approach that enables high-throughput discovery of polymorphic MHC class I-associated peptides (MIPs), which play a major role in allorecognition. On the basis of comprehensive analyses of the genomic landscape of MIPs eluted from B lymphoblasts of two MHC-identical siblings, we show that 0.5% of non-synonymous single nucleotide variations are represented in the MIP repertoire. The 34 polymorphic MIPs found in our subjects are encoded by bi-allelic loci with dominant and recessive alleles. Our analyses show that, at the population level, 12% of the MIP-coding exome is polymorphic. Our method provides fundamental insights into the relationship between the genomic self and the immune self and accelerates the discovery of polymorphic MIPs (also known as minor histocompatibility antigens).

Identification of small molecules that support human leukemia stem cell activity ex vivo.

Leukemic stem cells (LSCs) are considered a major cause of relapse in acute myeloid leukemia (AML). Defining pathways that control LSC self-renewal is crucial for a better understanding of underlying mechanisms and for the development of targeted therapies. However, currently available culture conditions do not prevent spontaneous differentiation of LSCs, which greatly limits the feasibility of cell-based assays. To overcome these constraints we conducted a high-throughput chemical screen and identified small molecules that inhibit differentiation and support LSC activity in vitro. Similar to reports with cord blood stem cells, several of these compounds suppressed the aryl-hydrocarbon receptor (AhR) pathway, which we show to be inactive in vivo and rapidly activated ex vivo in AML cells. We also identified a compound, UM729, that collaborates with AhR suppressors in preventing AML cell differentiation. Together, these findings provide newly defined culture conditions for improved ex vivo culture of primary human AML cells.

A key role for EZH2 and associated genes in mouse and human adult T-cell acute leukemia.

In this study, we show the high frequency of spontaneous γδ T-cell leukemia (T-ALL) occurrence in mice with biallelic deletion of enhancer of zeste homolog 2 (Ezh2). Tumor cells show little residual H3K27 trimethylation marks compared with controls. EZH2 is a component of the PRC2 Polycomb group protein complex, which is associated with DNA methyltransferases. Using next-generation sequencing, we identify alteration in gene expression levels of EZH2 and acquired mutations in PRC2-associated genes (DNMT3A and JARID2) in human adult T-ALL. Together, these studies document that deregulation of EZH2 and associated genes leads to the development of mouse, and likely human, T-ALL.

Microarray analysis of HIV resistant female sex workers reveal a gene expression signature pattern reminiscent of a lowered immune activation state.

To identify novel biomarkers for HIV-1 resistance, including pathways that may be critical in anti-HIV-1 vaccine design, we carried out a gene expression analysis on blood samples obtained from HIV-1 highly exposed seronegatives (HESN) from a commercial sex worker cohort in Nairobi and compared their profiles to HIV-1 negative controls. Whole blood samples were collected from 43 HIV-1 resistant sex workers and a similar number of controls. Total RNA was extracted and hybridized to the Affymetrix HUG 133 Plus 2.0 micro arrays (Affymetrix, Santa Clara CA). Output data was analysed through ArrayAssist software (Agilent, San Jose CA). More than 2,274 probe sets were differentially expressed in the HESN as compared to the control group (fold change ≥1.3; p value ≤0.0001, FDR <0.05). Unsupervised hierarchical clustering of the differentially expressed genes readily distinguished HESNs from controls. Pathway analysis through the KEGG signaling database revealed a majority of the impacted pathways (13 of 15, 87%) had genes that were significantly down regulated. The most down expressed pathways were glycolysis/gluconeogenesis, pentose phosphate, phosphatidyl inositol, natural killer cell cytotoxicity and T-cell receptor signaling. Ribosomal protein synthesis and tight junction genes were up regulated. We infer that the hallmark of HIV-1 resistance is down regulation of genes in key signaling pathways that HIV-1 depends on for infection.

Redefining interferon-producing killer dendritic cells as a novel intermediate in NK-cell differentiation.

The cell lineage origin of IFN-producing killer dendritic cells (IKDCs), which exhibit prominent antitumoral activity, has been subject to debate. Although IKDCs were first described as a cell type exhibiting both plasmacytoid DC and natural killer (NK) cell properties, the current view reflects that IKDCs merely represent activated NK cells expressing B220, which were thus renamed B220+ NK cells. Herein, we further investigate the lineage relation of B220+ NK cells with regard to other NK-cell subsets. We surprisingly find that, after adoptive transfer, B220- NK cells did not acquire B220 expression, even in the presence of potent activating stimuli. These findings strongly argue against the concept that B220+ NK cells are activated NK cells. Moreover, we unequivocally show that B220+ NK cells are highly proliferative and differentiate into mature NK cells after in vivo adoptive transfer. Additional phenotypic, functional, and transcriptional characterizations further define B220+ NK cells as immediate precursors to mature NK cells. The characterization of these novel attributes to B220+ NK cells will guide the identification of their ortholog in humans, contributing to the design of potent cancer immunotherapies.

Regulation of efflux pump expression and drug resistance by the transcription factors Mrr1, Upc2, and Cap1 in Candida albicans.

Constitutive overexpression of the Mdr1 efflux pump is an important mechanism of acquired drug resistance in the yeast Candida albicans. The zinc cluster transcription factor Mrr1 is a central regulator of MDR1 expression, but other transcription factors have also been implicated in MDR1 regulation. To better understand how MDR1-mediated drug resistance is achieved in this fungal pathogen, we studied the interdependence of Mrr1 and two other MDR1 regulators, Upc2 and Cap1, in the control of MDR1 expression. A mutated, constitutively active Mrr1 could upregulate MDR1 and confer drug resistance in the absence of Upc2 or Cap1. On the other hand, Upc2 containing a gain-of-function mutation only slightly activated the MDR1 promoter, and this activation depended on the presence of a functional MRR1 gene. In contrast, a C-terminally truncated, activated form of Cap1 could upregulate MDR1 in a partially Mrr1-independent fashion. The induction of MDR1 expression by toxic chemicals occurred independently of Upc2 but required the presence of Mrr1 and also partially depended on Cap1. Transcriptional profiling and in vivo DNA binding studies showed that a constitutively active Mrr1 binds to and upregulates most of its direct target genes in the presence or absence of Cap1. Therefore, Mrr1 and Cap1 cooperate in the environmental induction of MDR1 expression in wild-type C. albicans, but gain-of-function mutations in either of the two transcription factors can independently mediate efflux pump overexpression and drug resistance.

RNA-seq analysis of 2 closely related leukemia clones that differ in their self-renewal capacity.

The molecular mechanisms regulating self-renewal of leukemia stem cells remain poorly understood. Here we report the generation of 2 closely related leukemias created through the retroviral overexpression of Meis1 and Hoxa9. Despite their apparent common origin, these clonal leukemias exhibit enormous differences in stem cell frequency (from 1 in 1.4, FLA2; to 1 in 347, FLB1), suggesting that one of these leukemias undergoes nearly unlimited self-renewal divisions. Using next-generation RNA-sequencing, we characterized the transcriptomes of these phenotypically similar, but biologically distinct, leukemias, identifying hundreds of differentially expressed genes and a large number of structural differences (eg, alternative splicing and promoter usage). Focusing on ligand-receptor pairs, we observed high expression levels of Sdf1-Cxcr4; Jagged2-Notch2/1; Osm-Gp130; Scf-cKit; and Bmp15-Tgfb1/2. Interestingly, the integrin beta 2-like gene (Itgb2l) is both highly expressed and differentially expressed between our 2 leukemias (∼ 14-fold higher in FLA2 than FLB1). In addition, gene ontology analysis indicated G-protein-coupled receptor had a much higher proportion of differential expression (22%) compared with other classes (∼ 5%), suggesting a potential role regulating subtle changes in cellular behavior. These results provide the first comprehensive transcriptome analysis of a leukemia stem cell and document an unexpected level of transcriptome variation between phenotypically similar leukemic cells.